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Abstract:  The  paper  addresses  a  class  of 
problems  for  determining  the  acoustic 
interaction  of  time  harmonic  plane  waves 
impinging  upon  submerged  elastic  plates 
(fluid  backed  or  void  backed),  that  are  of 
infinite  in  extent.  The  plates  can  have 
repeated  equally  spaced  inclusions  (voids, 
solids,  ribs  etc.),  and  also  can  be 
constructed  from  orthotropic  and/or  isotropic 
layers  that  can  have  a  frequency  dependent 
loss  factor  for  the  modulus  data  (e.g.  E= 
E’(1+ir)E),  and  v=  v’(1+ir|v),  where  r|E  and  qv 
are  the  corresponding  loss  factors  for 
Young’s  modulus  and  Poisson’s  ratio  for 
example.  Since  the  plates  are  infinite  in 
extent,  the  issue  of  dealing  with  the  infinite 
domain  of  fluid  must  be  dealt  with  and  is 
treated  through  the  application  of  PML 
(Perfectly  Matched  Layers)  boundary 
conditions,  for  absorbing  the  reflected 
acoustic  waves  and  transmitted  acoustic 
waves  when  the  plates  are  fluid  backed.  The 
physical  quantities  of  interest  are  typically 
the  reflected  acoustic  waves  and  transmitted 
acoustic  waves,  however  detailed 
information  of  how  the  elastic  stress  waves 
propagate  through  and  around  the 
inclusions  are  often  of  interest  as  well. 
Keywords:  Fluid  Structure  Interaction,  Acoustic 
Transparency 

1.  Introduction 

The  paper  treats  a  general  class  of  problems 
as  encountered  in  structural  acoustics,  involving 
a  submerged  elastic  plate  that  is  subject  to  an 
incident  harmonic  plane  wave  of  the 

for ™  Pinc(x,y,z,t)  =  p0el{~xk*-yky-zK+<ot)  Eq(l) 
, where {kx  ky  kz  fare  the  acoustic  wave  number 
components  defining  the  direction  n  of 
propagation  (i.e.  n  =ikx  +  jkY  +  kkz)  and  co  is 

the  frequency  in  rad./sec  .  The  thrust  of  this 
paper  is  to  compute  the  acoustic  transparency 
|ptr/p0|  and  acoustic  reflection  |prt/p0|,  which  are 
defined  as  the  fractions  of  the  incident  pressure 
p0  that  are  transmitted  and  reflected  from  the 
plate  respectively.  Further,  we  are  interested  in 


the  case  where  the  plate  is  not  necessarily 
homogenous  in  the  plane  of  the  plate,  but  rather 
has  some  sort  of  repeated  anomaly,  such  as  the 
model  configuration  shown  in  Fig.  1  . 
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Figure  1..  Typical  Repeated  inclusion  Submerged 
Structure  Model  Subject  to  Incident  Plane  Wave 
(with  PML  Wave  Absorbing  Layers) 

The  paper  will  pass  through  a  sequence  of 
increasingly  difficult  problems  which  illustrate 
the  application  of  the  periodic  boundary 
condition,  starting  with: 

(a)  a  simple  free  field  block  of  fluid  (e.g.  Fig.  2) 

(b)  a  homogenous  submerged  elastic  plate  with 
no  inclusions,  (e.g.  Fig  1  case  but  with  the  softer 
surrounding  matrix  material  set  equal  to  the 
stiffer  inclusion  material) 

(c)  a  submerged  plate  with  periodic  inclusions 
for  the  problem  type  shown  in  Fig.  1  . 

Each  of  these  problems  has  a  specific  point 
to  illustrate,  namely:  problem  (a)  shows  the 
enforcement  of  the  periodic  boundary  condition 
in  the  cleanest  form  (no  fluid  structure 
interaction  to  deal  with),  and  employs  the  perfect 
(exact  in  this  case)  impedance  boundary 
condition  to  terminate  the  far  field;  problem  (b) 
illustrates  the  setup  for  fluid  structure  interaction 
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boundary  conditions,  absorbing  radiation 
boundaries,  and  the  issue  of  using  the  scattered 
pressure  formulation,  along  with  validation 
against  an  exact  solution,  and  finally  problem  (c) 
solves  the  problem  of  main  interest,  where 
inclusions  are  distributed  in  a  periodic  pattern. 

1.1  Past  Work 

The  motivation  for  this  work  stems  from  the 
need  to  evaluate  the  acoustic  transparency  of 
various  kinds  of  submerged  plates  and  shells, 
and  in  particular,  when  acoustic  sensors  are 
placed  behind  these  structures  and  their 
performance  depends  on  the  strength  of  the 
transmitted  signal  |ptr/p0|.The  closer  this  ratio  is 
to  1.0,  the  better  the  transparency  and  therefore 
the  better  the  performance.  There  is  a  related 
problem  of  target  strength  (but  will  not  be 
covered  herein),  where  it  is  desired  to  keep  the 
reflected  pressure  |prt/p0|  small,  usually  through 
energy  absorption  via  large  loss  factors  in  the 
material  elastic  constants.  For  either  the 
transparency  or  target  strength  application,  the 
methodology  addressed  herein  can  be  used  to 
attack  these  problems.  References  [  1  -4]  treat  the 
general  issue  of  acoustic  transparency  from  a 
materials  point  of  view  and  references  [5-7] 
employ  three  FEM  based  computer  codes 
(NASTRAN,  ATILA,  and  STARS3D 
respectively)  that  have  been  used  to  solve  for  the 
acoustic  response  of  submerged  periodic 
structures. 


2.  Governing  Equations 

The  governing  equations  for  the  total 
pressure  p  in  the  acoustic  domain  and 
displacement  vector,  u  ,• ,  in  the  solid  domain  (for 

eU0‘  time  harmonic  type  response)  are  given  by: 

c2(d2pldx 2  +d2pldy2  +d2p/dz2)  +  co2p  =  0 


c  2JL(^L)  +  c2JL(!“L 

cl  dx-ydx.’  s  dx-y  8x- 

J  l  l  l 


8u ■  7 

-^-L)  +  cozu  =0 

°Xj  J 


where  c,  cd  and  cs  are  the  acoustic  wave  speed, 
complex  solid  dilatational  wave  speed  and 
complex  solid  shear  wave  speed  respectively. 
The  live  index  j  takes  on  values  j=l,2,3  and 
repeated  indices  sum  over  1  to  3.  We  write  the 
governing  equations  in  this  form,  in  order  to:  (a) 
illustrate  the  similarity  of  the  form  of  the  solid 
and  acoustic  media,  and  (b)  illustrate  how  wave 
speeds  appear  as  primary  coefficients  in  the 


harmonic  field  equations.  The  concept  of 
acoustic  transparency  in  the  simple  case  of  say  a 
normal  incident  plane  wave,  Eq(l),  impinging  on 
a  homogenous  elastic  plate  without  loss,  is 
treated  in  Ref  [8].  The  normal  incident  wave  will 
completely  pass  through  the  plate  if  the 
impedance  ratio  Z=pscd/pc  =1.0,  where 
p  and  ps  are  the  fluid  and  solid  plate  mass 
density  respectively.  Therefore  when  selecting 
material  parameters  for  a  plate,  where  acoustic 
transparency  is  desired,  one  usually  strives  to 
make  the  plate  materials  have  impedance  ratios 
in  the  neighborhood  of  1.0  .  In  the  case  of  non 
normal  incidence,  and/or  the  case  where  the 
matrix  material  of  the  plate  has  elastic  inclusions 
imbedded  within  the  plate  and/or  has  nonzero 
loss  factors,  other  kinds  of  waves  are  present  in 
the  dynamic  response,  and  then  having  Z  » 1 .0  to 
achieve  transparency  is  only  used  as  a  guideline. 


2.1  Elastic  Material  Constants 

The  complex  wave  speeds  can  be  expressed  in 
terms  of  primed  real  wave  speeds  and  loss 
factors  in  the  form: 


c\  =(0)2(!  +  /^) 


Cs  =(c's  >2(1  +  ‘V 


Eqs(2a) 


Where  r)j  and  r|s  are  the  dilatational  and  shear 
loss  factors.  Complex  wave  speeds  can  also  be 
expressed  in  terms  of  complex  E  modulus  and 
complex  Poisson’s  ratio  v  in  the  form: 


c2  (1  ~  v)E  c2 

d  Ps  (1  -  2v)(l  +  v) 


0.5  E 
Ps(  1  +  v) 


Eqs(2b) 


Say  we  are  given  the  two  real  primed  wave 
speeds  and  two  associated  loss  factors  r|d  and  r|s 
(which  are  typically  different  for  rubber  like 
materials).  Upon  equating  the  corresponding  like 
complex  wave  speeds  of  Eqs(2a)  Eqs(2b),  we 
can  easily  solve  for  complex  E  and  complex  v  in 
terms  of  c'j  ,  c'  ,  iqd  and  r|s  .  In  COMSOL,  the 

isotropic  elastic  data  is  entered  via  E  and 
v  ,  where  one  can  enter  the  complex  E  and  v  in 
the  form  E=E’  +  i*  E”  and  v=v’  +  i*  v”.  The 


COMSOL  program  has  a  feature  in  the  acoustics 
module  that  permits  the  user  to  work  with  the 
scattered  (reflected)  pressure  formulation,  ps, 
rather  than  the  total  p  formulation.  Flowever  the 
formulation  used  here  still  uses  the  structure  of 


the  total  pressure  equations,  where  the  loading  of 
the  incident  field  enters  via  the  boundary  loading 
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at  the  incident  side  fluid  structure  interface.  We 
start  by  substituting  p  =  pinc+  ps  into  the  original 
time  harmonic  pressure  pde  to  get: 

c“(0“ p$  ! 8x ^  +d~ps  / dy~  +  d^ps  ldz~)  +  arps  =0 

and  since  p;nc  is  a  solution  to  the  original 
homogenous  pde,  it  drops  out  after  the 
substitution,  leaving  the  new  form  of  the  pde  in 
terms  of  ps  which  has  the  same  appearance  as 
the  original. 

2.2  Fluid  Solid  Interface  Boundary  Conditions 

The  pinc  term  enters  the  problem  via  the  fluid 
structure  interface.  On  the  solid  interface,  (pjnc+ 
ps)=  the  total  pressure  is  applied,  and  similarly  at 
the  fluid  interface,  both  the  interface  plate 
acceleration  and  incident  fluid  acceleration  are 
loaded.  For  example: 

On  incident  side  interface: 
inward  normal  acceleration  fluid  loading : 
nx_acpn*u_tt_acpn+ny_acpn*v_tt_acpn+p_i_ac 
pr*i*(nx_acpn*kx+ny_acpn*ky)/rhow 
normal  pressure  loading 
-p*ny_acpn-p*nx_acpn-p_i_acpr*ny_acpn- 
p_i_acpr*nx_acpn 

On  transmitted  side  interface: 
inward  normal  acceleration  fluid  loading : 
nx_acpn*u_tt_acpn+ny_acpn*v_tt_acpn 
normal  pressure  loading 
-p  *ny_acpn-p  *nx_acpn 

Where  p_i_acpr  =  Eq(  1 )  without  e1®';  rhow  =  p 
po=1.0,  kx=-k0sin(9),  ky=k0cos(9),  k z  0  k.^m/c  . 

2.3  Radiation  Boundary  Conditions 

The  mesh  termination  at  the  top  of  the  incident 
side  fluid  and  bottom  of  the  transmitted  side 
fluid  must  include  some  sort  of  radiation 
absorbing  boundary  condition.  Three  types  were 
considered: 

a)  COMSOL ’s  built  in  plane  wave  radiation 
absorber  (don’t  need  acoustics  module') 

b)  COMSOL’s  built  in  PML  (Perfectly  Matched 
Layer),  as  indicated  by  the  extra  PML  zones 
shown  in  Fig.l  above  (  need  acoustics  module). 

c)  User  defined  Impedance  (Z=ps/vn)  radiation 
condition  ,  where  for  plane  waves  impinging  at 


angle  9  (see  Fig.  1),  on  one  or  more  completely 
submerged  homogenous  elastic  layers  and  using 
Snell’s  law,  ref.[9],  we  have  Z=pfc/cos(9)  as  an 
exact  radiation  boundary  condition  at  both 
incident  side  mesh  termination  and  transmitted 
side  termination.  This  same  impedance  condition 
can  be  used  for  the  non-homogenous  plate  (e.g. 
with  inclusions),  however  in  this  case,  the 
condition  is  approximate  and  the  degree  of 
success  depends  on  how  close  the  solution  field 
is  to  a  plane  wave  structure,  at  the  fluid  domain 
truncation  cut  (don ’t  need  acoustics  module). 

2.4  Periodic  Boundary  Conditions 

Here  we  are  concerned  about  the  boundary 
condition  at  the  left  cut  vertical  faces  and  right 
cut  vertical  faces  (line  x=0  and  line  x=L  in  Fig. 
1.  The  response  at  the  left  cut  say,  picut,  and 
response  at  the  right  cut  prcut,  is  not  known  in 
advance,  however  we  do  know  a  relationship 
between  them,  namely: 

Prcut  =  Picute'kxL ,  therefore  pkut  and  prcut  are  not 

independent  unknowns.  This  type  of  condition  is 
used  in  refs. [5-7]  for  example,  and  will  not 
elaborate  on  them  further.  Applying  this 
condition  in  COMSOL  is  awkward  and  not 
intuitive.  There  are  two  issues:  a)  when  applying 
the  periodic  boundary  involving  a  complex 
multiplier  condition,  the  user  applies  variable  p 

on  the  left  and  p  e~l  ,L  on  the  right  (the 

conjugate  of  e,kxL  )  and  b)  when  using  the  solve 
parameters  advanced  settings,  one  must 
(according  to  COMSOL  staff),  check  the  box 
engaging  “Use  Hermitian  transpose  of  constraint 
matrix). 

The  treatment  for  the  application  of  the  periodic 
boundary  condition  for  the  solid  follows  exactly 
along  similar  lines,  (where  upon  defining 
displacement  components  as  u=Ui  and  v=u2)  we 
have  urcut=ulcuteik*L  and  vrcut  =  vlcuteik*L  with 
{  ,  view  ,  urcut ,  v rcut }  as  the  left  and  right  cut 

displacement  component  values,  analogous  to 
the  pressure  left  and  right  cut  values.  The 
displacements  are  enforced  similar  to  enforcing 
p,  but  instead  in  terms  of  COMSOL  variables 
{u,v}.  Thus  {u,v}  are  applied  at  the  left  cut,  and 

{u  e,k,:L  ,  velk'-L  }are  applied  at  the  right  cut. 

2.5  Mesh  sizing 

For  harmonic  steady  state  problems,  the  mesh 
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size  is  set  according  to  the  shortest  wave  length 
expected  during  the  event.  For  example,  if  New  is 
the  number  of  elements/wave  length  required  for 
accurate  modeling  and  Cmjn  is  the  slowest  wave 
speed,  and  fmax  the  largest  frequency  experienced 
in  a  frequency  sweep,  then  the  mesh  can  be 
sized  with  Anljn  —  Cmjn  /  (New  fmax  ),  ( e . g.  New  6 
for  quadratic  element  shape  functions  and  New  = 
10  for  linear  element  shape  functions).  When 
rubber  like  materials  are  employed,  the  shear 
wave  speed,  c's ,  is  typically  the  smallest  and 
governs  the  mesh  size  needed  in  the  solid. 

2.6  Post  Processing  Pressure  Fields 

When  the  scattered  formulation  is  employed  as 
used  in  this  paper,  while  plotting  the  pressure  on 
the  incident  side,  the  basic  COMSOL  variable 
p  t  acpr  will  actually  represent  the  scattered 
component,  ps,  even  though  it  is  labeled  as  total 
in  the  post  processing  output  list.  The  actual  total 
can  be  obtained  (if  needed)  by  simply  adding 
back  the  incident  pressure,  Eq(l),  (without  elmt). 
The  pressure  on  the  transmitted  side  is  also  post- 
processed  via  COMSOL  via  the  variable 
p_t_acpr. 

When  the  plate  is  homogenous,  the  phase  of  the 
reflected  and  transmitted  pressure  vary  with 
space,  however  the  magnitudes  of  these  fields 
are  constant  over  the  domain.  Therefore,  one  can 
sample  |ptI/p0|  ,  or  |pri/p0|  anywhere  in  the  field  to 
get  a  representative  value  of  the  acoustic 
transparency  and  reflectivity.  However  in  the 
case  of  the  presence  of  a  repeated  inclusion,  like 
in  the  Fig.l  model,  one  needs  a  strategy  for 
computing  the  transmitted  and  reflected  pressure. 
One  simple  measure  would  be  to  compute 
integrated  pressure  at  the  far  field  mesh 
boundary  (but  not  inside  the  PML  zone  if 
present),  thus  getting: 


where  integration  along  x  is  at  y=constant  (the 
fluid  mesh  boundary). 

Next,  we  consider  an  alternate  method  to 
process  a  representative  transparency  pressure  in 
a  variable  spatial  field,  namely  by  computing  the 
power  flow  across  the  y=constant  cut  boundary, 
and  then  convert  the  power  into  an  equivalent 
pressure  we  shall  call  p  =  "pseudo  pressure ”. 


The  power  flow  can  be  computed  by  integrating 
the  work  done  over  one  time  cycle,  (e.g.  as  in 
Ref.[10]),  and  then  integrating  that  power/area 
result  over  the  top  (or  bottom  )  boundary  cut 
surface  of  the  finite  element  mesh,  getting 

L 


n=  - 


J^Re(pvn*)& 


, where  v*n  denotes  the 


o 

complex  conjugate  of  the  velocity  normal  to  the 
surface  .  The  power  of  the  incident  wave  is  given 


by  Ylinc  =(\l2)plLcos{0)l(cp f )  and  is  used  to 


form  the  normalized  power  ratio  for  transmitted 
and  reflected  fields,  which  we  defined  as 
n„-  —  n !r  n and  nr/  —  n,y/n7-HC  .  Eqs(4) 


These  expressions  can  be  used  to  compute  the 
pseudo  pressure  ratio  \  ptr  |  /  pa  ,  by  equating  the 


normalized  power  Htr ,  to  the  normalized  power 
in  a  pseudo  plane  wave  (i.e. 

If,,.  =  (1  /  2)piL  cos(0)  i(cpf )  IYlinc ),  and  solving 

for  the  transmitted  pseudo  pressure  (and 

similarly  the  reflected  value),  we  obtain 

\Ptr\l  Po=^\-  \Prf  \/Po  •  Eqs(5) 

When  the  plate  is  homogenous,  the  reflected  and 
transmitted  waves  are  plane  with  constant 
amplitudes,  then  Eqs(3)  and  Eqs(5)  yield  the 
same  value  for  normalized  pressure. 


The  power  relations  of  Eqs(4)  have  an  alternate 
use,  namely  for  checking  the  consistency  of  the 
FEM  solutions.  Power  levels  can’t  be  created 
greater  than  the  input  normalized  incident  power 
level  of  n,„c=l,  hence  the  following  inequality 
must  be  met: 

n„.  +  nr/<nj„c=l.O  .  Eq(6) 

The  strict  equality  sign  is  used  when  all  the  plate 
loss  factors  are  zero.  When  the  plate  loss  factors 
are  not  zero,  the  inequality  is  used  due  to 
material  energy  dissipation,  where  Eq(5)  then 
infers  that  the  transmitted  plus  reflected 
normalized  power  is  less  than  the  normalized 
free  field  incident  power  (i.e.  1.0).  When 
computing  power  with  Eq(4)  at  the  fluid-PML 
interface,  the  vn  quantities  must  be  computed 
with  one  sided  shape  function  data  (in  the  regular 
acoustic  side),  by  using  the  </ow/i(vy_acpr) 
option  for  the  reflected  side  fluid  and 
up(v y  aepr)  option  for  the  transmitted  side  fluid. 
Otherwise  accurate  surface  velocities  will  not  be 
obtained. 
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3.  Applications 

Here  we  give  examples  of  the  three  problem 
types  (a,b,c)  outlined  in  the  introduction  section. 

3.1  Free  Field  Propagation  Through  a  Simple 
Fluid  Block 

The  purpose  of  this  demonstration  is  to  illustrate 
the  enforcement  of  the  periodic  boundary 
conditions,  with  out  the  additional  complications 
of  PML  absorbers,  and  fluid  structure  interaction 
boundary  conditions.  The  model  is  shown  in  Fig. 
2,  where  the  bottom  is  driven  with  Eq(l) 
(without  elmt  since  it  is  suppressed  throughout  the 
harmonic  solution),  with  kz=0  and  k0=co/c,  where 
positive  kx  and  ky  values  send  the  wave  in  the 
direction  shown  in  Fig.  2. 

In  this  special  example  we  work  directly  with 
total  pressure  p,  since  there  is  no  issue  of 
reflections  in  the  solution.  Further,  the  model  is 
terminated  with  the  perfect  impedance  absorber 
Z=pfc/cos(0),  so  there  is  no  need  for  PML 
absorbers  here. 

The  left  and  right  boundary  conditions  are 
treated  by  two  methods: 

(i)  a  direct  Dirichlet  type  b.c.  where  the  left 
and  right  cuts  are  directly  driven  with  Eq(l) 
at  x=0  for  the  left  cut  and  with  Eq(l)  with 
x=L  at  the  right  cut  (without  elmt). 

(ii)  a  periodic  boundary  condition  was 
employed  to  the  left  and  right  vertical  cuts, 
and  since  there  is  no  characteristic  period 
length,  L  can  be  set  as  anything,  thus  a 
nominal  L=1.0”  was  selected.  The  left  and 
right  cuts  were  driven  by  the  method 
explained  in  section  2.4,  with  strict  attention 
to  the  issues  (a)  and  (b)  mentioned  therein. 

A  l”x  6”  block  of  fluid  was  modeled,  and  driven 
at  a  frequency  f=10000  Hz,  at  an  angle  0=60°. 
The  water  wave  speed  was  c=60000  in/sec  and 
the  mass  density  pf=.000096  sec-lbf/in4  and  the 
mesh  size  was  set  via  section  2.5  . 

The  solutions  for  the  (i)  and  (ii)  methods  of 
driving  the  vertical  faces  were  the  same,  thus 
validating  the  procedure  for  enforcing  the 
periodic  boundary  condition.  The  solution  for  the 
(ii)  case  is  shown  in  Fig. 3  and  perfectly  checks 
out  against  the  exact  solution  (i.e.  the  |p| 
magnitude  should  be  1.0  in  the  entire  field,  and 
the  direction  of  wave  propagation  should  be 


along  straight  lines  that  are  at  angles  9  =60°  to 
the  y  axis,  where  it  is  noted  that  these  lines  are 


Z=pC/COS(0) 


2  D  model 


(inf.  out  of  plane) 
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Figure  2.  FEM  model  Zone  (dashed)  for  Free  Field 
Propagating  Plane  Wave  Traveling  0  Degrees  to 
the  y  Axis. 


perpendicular  to  the  lines  of  constant  phase  wave 
fronts.  This  clean  simple  demo  validates  the 
procedure  for  enforcing  the  periodic  boundary 
condition. 


Mat  l  BO  1  DO 


MM  140  0 


Figure  3.  FEM  Pressure  Magnitude  and  Phase  for 
Free  Field  Propagating  Plane  Wave  Traveling 
9  =60  0  to  the  y  Axis. 
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3.2  Homogenous  Submerged  Elastic  Plate 
with  No  Inclusions  (with  and  without  loss) 

This  example  is  selected  in  order  to  validate  the 
solution  in  the  more  complicated  case,  where  in 
addition  to  the  periodic  boundary  condition,  the 
issue  of  dealing  with  an  elastic  plate  (with  and 
without  loss)  is  treated.  The  model  setup  is 
shown  in  Fig.  1,  except  the  softer  surrounding 
matrix  material  set  equal  to  the  stiffer  inclusion 
material),  thus  ending  up  with  a  single  (all 
orange)  material  plate. 

Since  the  elastic  layer  is  now  simply  flat  and 
homogenous,  we  can  directly  compare  the  FEM 
solutions  to  an  exact  solution  of  the  same 
problem  taken  from  Ref.  [9]  .  The  model  consist 
of  a  5.5”  wide  x  6”  long  front  and  back  side 
fluid,  truncated  with  an  optional  PML  zone  of 
5.5”x  4”  .  The  water  properties  are  the  same  as 
the  previous  problem  and  the  drive  frequency  is 
f=l  1,000  Hz.  The  solid  properties  have  a 
dilatational  wave  speed  magnitude  in  the  general 
vicinity  of  water,  but  still  has  a  substantial  shear 
modulus,  hence  a  substantial  shear  wave  speed. 
More  specifically,  the  for  the  plate  we  have: 
c’d  =  80000  in/sec,  c'  =20000  in/sec, 

qd=0.01,  qs  =  0.10,  and  ps=. 00011  sec-lbf/in4, 
where  upon  using  the  relations  explained  in 
section  2.1,  corresponds  to  direct  COMSOL 
domain  input  E=  129,090  +i*126250  and 
v=0. 46665  -i  *0.031999  .  This  material  has  an 
impedance  mismatch  ratio  of  Z  =  pscd!  pc= 1.52 
,  which  is  well  away  from  the  ideal  1.0  and 
therefore  it  is  expected  that  unwanted  reflections 
appear  as  will  be  shown  later  in  Fig. 4  .  When 
entering  this  modulus  and  Poisson’s  ratio  data, 
we  skip  COMSOL’s  built  in  loss  mechanism 
input  which  assigns  the  same  loss  (unwanted 
here)  to  both  dilatation  and  shear.  Instead,  we 
enter  the  data  directly  as  MATLAB  like  syntax 
entries.  The  frequency  is  held  constant  and  swept 
over  a  range  of  angles  of  incidence  0<9<60°. 
The  same  problem  was  solved  for  three  radiation 
boundary  absorbers  described  in  section  2.3  and 
FEM  solutions  are  compared  to  the  exact 
solution  (for  transmitted  pressure,  reflected 
pressure  and  Power).  The  simpler  Eqs(3)  were 
used  to  process  pressure  across  the  mesh 
truncating  fluids  faces,  and  Eqs(4)  for  the  power. 
The  comparisons  in  Fig  4.  are  excellent  for  both 
the  PML  and  the  Z=pfc/cos(0)  absorbers.  Note: 
in  all  plots,  “MPL”  labels  imply  PML  absorbers. 
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Figure  4.  Exact  vs.  FEM  Solution  for  Three 
Different  Radiation  Boundary  Conditions  (With 
Loss ,  f=ll  KHz). 


The  same  problem  was  solved  for  no  loss  in  the 
solid  (setting  the  imaginary  parts  of  E  and 
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Figure  5.  Exact  vs.  FEM  Solution  for  Three 
Different  Radiation  Boundary  Conditions  (With 
No  Loss,  f=ll  KHz). 
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v  equal  to  zero),  where  the  results  are  shown  in 
Fig.  5  .  We  note  that  since  the  loss  is  zero,  the 
total  power,  via  Eq(6),  should  be  conserved  (i.e. 
be  equal  1.0  all  angles).  However  as  in  the 
previous  plot,  the  power  computed  with  the 
COMSOL  plane  wave  absorber,  drifts  off  the 
exact  solution  of  1 .0  for  30°  <0<  60°  . 

3.3  Non-homogenous  Submerged  Elastic  Plate 
With  Periodic  Inclusions  (with  and  without 
loss) 

The  model  is  the  same  as  the  Fig.  1  sketch  (drawn 
to  scale),  except  the  matrix  material  is  a  rubber 
like  material,  and  a  centered  l”x  1”  inclusion 
made  of  the  same  material  as  the  previous 
homogenous  plate  example.  The  matrix  material 
(green),  is  given  as:  c'^=  65000.  in/sec, 

c'  =10000.,  r|d  =0.03,  r|s  =  0.30  ,  p  =.000105 

sec-lbf/in4,  which  upon  using  the  relations 
explained  in  section  2.1,  corresponds  to  direct 
COMSOL  domain  input  E=  31265  +i*9302.8 
and  v=0.4878  -i  *0.0033502.  This  matrix 
material  has  an  impedance  mismatch  ratio  of 
Z matrix  =  Pscd!  P°=  1-18  ,  which  is  near  the  ideal 
1.0,  however  the  stiffer  inclusions 
(Z infusion— 1.52  )  are  expected  to  somewhat 
degrade  the  overall  transparency  performance. 


3.3.1  Baseline  case:  The  resulting  solution  is 
shown  in  Fig  6,  where  note  that  results  using 


«  ANGLE  OF  INCIDENCE-;.;  FREG=1 1.00  KHZ 

Figure  6.  FEM  Solution  for  Three  Different 
Radiation  Boundary  Cond.  (With  Loss,  f=ll  KHz) 


the  PML  and  Z  impedance  absorbers  are  in  good 
agreement,  however  the  COMSOL  plane  wave 
absorber  solution  drifts  away  for  30°  <0<  60°. 


Next  we  consider  the  same  solution,  except  the 
loss  terms  are  turned  off  for  both  the  matrix 


material  and  the  inclusion.  The  resulting  solution 
is  shown  in  Fig.  7,  where  this  no  loss  case  is 
included  in  order  to  illustrate  how  the  Eq(6) 
power  inequality  checks  the  solution  accuracy. 
Thus  for  the  COMSOL  built  in  plane  wave 
absorber,  the  power  sum  does  not  always  add  up 
to  1.0,  whereas  with  both  the  PML  and  the  Z=pf 
c/cos(0)  absorbers,  it  does  nicely  add  up  to  1.0 
across  the  entire  angle  sweep,  even  when  the 
solution  encounters  the  spikes  in  the  response. 
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Figure  7.  FEM  Solution  for  Three  Different 
Radiation  Boundary  Cond.  (No  Loss,  f=ll  KHz). 


The  color  contour  plots  in  Fig.  8,  illustrate  how 
the  transmitted  and  reflected  pressures  vary  with 
the  spatial  x,y  field  ,and  how  phase  wave  fronts 
exist  in  the  transmitted  side  and  reflected  side 
within  the  fluid  domain  at  0=  60°. 


3.3.2  Variation  of  Baseline  Case  for  three 
Inclusion  Stiffnesses: 

In  this  demonstration,  we  illustrate  how  the 
transmitted  pressure  is  affected  by  the  stiffness 
of  the  inclusion.  The  stiffnesses  are  changed  via 
the  solid  dilatational  wave  speeds.  The  without 
inclusion  case  is  run  first  for  comparison 
purposes,  where  the  entire  plate  is  made  of  the 
matrix  material  (red  curve  solution).  This  is  used 
as  a  yardstick  reference,  for  measuring  the 
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Figure  8.  FEM  Amplitude  and  Phase  Solution  for 
PML  Radiation  Boundary  Conditions  (With  No 
Loss,  f=ll  KHz,  at  0=  60°). 

without  inclusion  solution  against  the  with 
inclusion  solution.  The  blue  curve  represents  the 
with  inclusion  base  line  case  for  the  case,  c' ^  = 
80000  in/sec  (Fig. 6  repeat),  and  the  green  curve 
for  the  stiffest  inclusion  with  c'd=  160,000 

in/sec  (keeping  all  other  material  parameters  the 
same  as  the  baseline  case).  The  sequence  of 
inclusion  stiffnesses  for  the  (Fig. 9  red  blue  green 
curves)  are  Z  ={1.0,  1.18,  3.05}.  As  expected, 
Fig. 9,  illustrates  that  the  stiffer  the  inclusion 
material,  the  worst  the  acoustic  transparency. 


INCLUSION  STIFFNESS  VARIATIONS  WITH  LOSS 


«  ANGLE  OF  INCIDENCE^;  FREQ=11jM  KHZ 

Figure  9.  Variation  of  Baseline  Case  for  Three 
Inclusion  Stiffnesses  (With  Loss,  f=ll  KHz). 


3.3.3  Variation  of  Baseline  Case  for  Two 
Fluid  Domain  Lengths  and  Two  Pressure  Post 
Processing  Techniques: 

In  this  example,  we  examine  the  effect  of  the 
fluid  length  on  the  accuracy  of  the  FEM  solution 
for  the  baseline  case  without  loss  and  also  the 
effect  of  average  vs.  pseudo  pressure  post¬ 
processing  as  shown  in  Fig.  10.  We  purposely 
zero  out  the  loss  portion  of  the  material 
constants,  since  there  is  a  richer  profile  of 
response  vs.  0  in  order  to  better  examine  the 

PROC  ESS  P  VIA  AVERAGE  VS  VIA  POW  ER  NO  LOSS 
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Figure  10.  Variation  of  Baseline  Case  for  Two 
Fluid  Lengths  and  Two  Post-processing  Methods 
(No  Loss,  PML  absorber  at  f=ll  KHz  ). 

pressure  post-processing  techniques  under  more 
severe  spatial  distribution  conditions.  The  curves 
labeled  MPL  LONG  refer  to  the  6“  length  front 
and  back  side  fluid  domains,  and  the  curves 
labeled  MPL  SHORT  refer  to  a  0.60”  length 
front  and  back  side  acoustic  domain.  The 
detailed  pressure  contours  (not  shown  here), 
indicated  that  the  shorter  fluid  domain  cuts  are 
now  in  the  near  field,  where  the  solution  pressure 
amplitude  vary  with  x  more  substantially  than 
the  long  model,  and  thereby  provide  a  good 
opportunity  to  compare  average  pressure,  Eq(3), 
vs.  pseudo  power  based  pressure  ,Eq(5),  post¬ 
processing.  The  results  in  Fig.  10,  illustrate  the 
pseudo  pressure  of  the  short  model  and  the  long 
model  produce  the  same  measure  of  transmitted 
pressure  and  reflected  pressure,  implying  one  can 
employ  a  shorter  model  (and  hence  less  DOF  in 
the  event  core  memory  is  an  issue  for  higher 
frequency  models).  However,  the  average 
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pressure  from  Eq(3),  (dashed  red)  of  the  shorter 
model  reflected  pressure  does  not  completely 
track  the  corresponding  long  model  pseudo 
pressure,  missing  at  peak  0=24°  and  valleys 
0=32°,  37°  .  This  illustrates  that  the  pseudo 
pressure  post-processing  is  a  more  general 
procedure  for  evaluating  the  transparency  and 
reflectivity  of  panel  type  materials  with 
inclusions,  and  that  it  produces  results  for  short 
length  models  that  are  consistent  with  the  longer 
(larger  DOF)  models. 

Finally  we  consider  the  same  long-short  model 
comparison,  except  now  the  loss  factors  are 
turned  on  for  both  materials,  resulting  in  the 
solution  shown  in  Fig.  11.  Unlike  the  previous 
Fig.  10  no  loss  case,  we  have  essentially 
complete  agreement,  across  the  entire  0  sweep, 
between  the  LONG,  SHORT  model  solutions, 
and  also  agreement  between  the  method  of  post¬ 
processing.  The  one  exception  is  a  relatively 
small  reflected  pressure  blip  (red  dashed  curve), 
near  normal  incidence  using  Eq(3)  pressure 
averaging. 


Figure  11.  Variation  of  Baseline  Case  for  Two 
Fluid  Lengths  and  Two  Post-processing 
Techniques  (With  Loss,  PML  Absorber  at  f=ll 
KHz). 

4.  Conclusions 

The  results  in  this  paper  illustrate  how 
COMSOL  could  be  used  to  solve  for  the 
transmitted  and  reflected  pressure  in  totally 
submerged  plates,  having  periodically  spaced 


inclusions.  Although  we  have  emphasized  the 
issue  of  transparency  with  simple  inclusions,  the 
methodology  could  be  applied  to  ribbed  plates, 
where  the  periodically  spaced  ribs  would  play 
the  same  role  as  the  different  material  inclusions 
considered  herein.  Also  we  have  shown  another 
method  for  post-processing  a  variable  spatial 
distributed  pressure  field,  by  computing  a  pseudo 
pressure,  derived  from  the  amount  of  power 
flowing  across  the  cuts  in  the  fluid  domain  outer 
boundaries.  The  power  flow  calculations  also 
provide  a  measure  ,via  Eq(6),  of  the  accuracy  of 
the  FEM  model  solution  .  Violations  of  Eq(6), 
provide  a  red  flag  that  something  is  wrong  with 
the  model  (e.g.  coarse  mesh,  bad  absorbing 
boundary,  some  unknown  user  error). 
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PAPER  OVERVIEW 


*  The  paper  addresses  a  class  of  problems  for  determining  the 
acoustic  interaction  of  eicot  time  harmonic  plane  waves  impinging 
upon  submerged  elastic  plates  (fluid  backed  or  void  backed),  that 
are  of  infinite  in  extent. 

*  The  plates  can  have  repeated  equally  spaced  inclusions  (voids, 
solids,  ribs  etc.),  and  also  can  be  constructed  from  orthotropic 
and/or  isotropic  layers 

9  Materials  can  have  a  frequency  dependent  loss  factor  for  the 
modulus  data  (e.g.  E=  E’(1+ir|E ),  and  u=  o’(1  +ir|u),  where  r|Eand  r|0 
are  the  corresponding  loss  factors  for  Young’s  modulus  and 
Poisson’s  ratio  for  example. 

*  THE  ACOUSTIC  TRANSPARENCY :  I.E.  INSERTION 
LOSS=MAG(Ptr/P|NC)  ,  is  the  main  physical  quantity  of  interest 

9  Transparency  is  of  great  interest  when  sonar  devices  operate 
behind  protective  dome  enclosures  (“  acoustic  windows”) 
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DEFINITION  OF  ACOUSTIC  WINDOW  KEY  PARAMETERS 


PAPER  OVERVIEW  cont. 


*  Since  the  plates  (and  surrounding  fluid)  are  infinite  in  extent,  the 
issue  of  dealing  with  the  infinite  domain  of  fluid  must  be  dealt  with 
and  is  treated  through  the  application  of  wave  absorbing 
boundary  conditions  for  absorbing  the  reflected  and  transmitted 
acoustic  waves 

*  Three  absorbing  boundary  conditions  are  examined 

(1)  COMSOL’s  PML  (Perfectly  Matched  Layers)  boundary 
conditions 

(2)  Impedance  boundary  condition 

(3)  COMSOL’s  plane  wave  absorbing  boundary  condition 

9  Only  a  unit  cell  of  the  repeating  structure,  including  a  small  part  of 
the  surrounding  fluid  domain  has  to  be  modeled,  by  applying  the 
the  Bloch-Floquet  theorem  based  boundary  condition  in  the 
direction  of  the  repeat  pattern 


FIG.  1  _Typical  Repeated  inclusion  Submerged  Structure  Model  Subject 
to  Incident  Plane  Wave  (with  PML  Wave  Absorbing  Layers) 


Past  Work  on  repeated  discontinuities  using,  periodic 
structures 

•  Ref-1  “Plates  with  regular  stiffening  in  acoustic  media: 
Vibration  and  radiation”,  Mead  D.J.,  1990,  J.  Acoustic  Society  of 
America,  88  (1),  July  1990. 

This  important  paper  treats  plates  and  cylindrical  shells  with 
regular  and  identical  stiffening  that  constitute  “spatially  periodic 
structures”.Specially  convenient  methods  of  vibration  analysis 
are  considered  for  these  cases,  some  of  which  are  suitable  for 
the  inclusion  of  the  effects  of  fluid  loading  from  adjacent 
acoustic  media. 

The  paper  outlines  the  nature  of  free  wave  motion  in  periodic 
structures  stiffened  either  in  one  or  two  orthogonal  directions. 

In  our  application,  we  could  view  the  stiffened  sections  say  due 
to  ribs  for  example,  as  unstiffened  -or-weakened  stiffness 
sectors,  like  the  failed  region  LF  we  illustrated  earlier. 


Past  Work  on  repeated  discontinuity  using  periodic  structures  cont. 

•  Ref-2  “Steady  state  solutions  to  dynamically  loaded  periodic 
structures”, Kalinowski  A.J.  1979,  NASA  Conference  Publication  2131, 
Eight  NASTRAN  User’s  Colloquium. 

This  was  an  early  application  of  periodic  structures  applied  to  acoustic 
scattering  and  transmission  of  acoustic  structure  interaction  problems 
having  one  way  periodicity  based  early  Mead  papers  predating  his 
general  Ref-29  paper.  DMAP  matrix  operations  in  the  NASTRAN  FEM 
code  were  used  for  performing  matrix  manipulations  needed  to  extract 
acoustic  solutions.  The  beauty  of  this  approach  is  that  one  can  interact 
with  the  inner  workings  of  the  NASTRAN  FEM  code  without  having  to 
have  access  to  the  source  code.  This  is  a  rare  feature  not  readily  found 
in  today’s  large  scale  FEM  codes. 
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Past  Work  on  repeated  discontinuity  using  periodic  structures  cont. 

•Ref-3  “Analysis  of  the  scattering  of  a  plane  acoustic  wave  by  a 
periodic  elastic  structure  using  the  finite  element  method:  Application 
to  compliant  tube  gratings”,  Hen n ion  A.C.,  Bossut  R.  and  Decarpigny 

J.N.,  1990,  J.  Acoustic  Society  of  America,  87  (5),  May  1990. 

A  mathematical  model  has  been  developed  to  analyze  the  scattering  of 
plane  acoustic  waves  from  an  infinite,  uniform,  plane  grating  of 
compliant  tubes.  It  relies  upon  the  FEM  method  and  uses  the  ATI  LA 
FEM  code  (see  Ref  30).  Only  a  unit  cell  of  the  repeating  structure, 
including  a  small  part  of  the  surrounding  fluid  domain  has  to  be 
modeled,  by  relying  on  the  Bloch-Floquet  theorem,  and  the  effects  of 
the  remaining  fluid  domain  are  accounted  for  by  matching  the  pressure 
field  in  the  finite  element  model  with  the  simple  plane  wave  expansions 
of  the  on  going  and  outgoing  waves. 
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Past  Work  on  repeated  discontinuity  using  periodic  structures 
cont 

•  Ref-4  “Analysis  of  the  scattering  of  a  plane  acoustic  wave  by  a 
doubly  periodic  structure  using  the  finite  element  method: 
Application  to  Alberich  anechoic  coatings”,  Hlady-Hennion  A.C., 
andDecarpigny  J.N.,  1991, Dec,  J.  Acoustic  Society  of  America, 
90(6). 

Ref-5”Computation  of  Acoustic  Transmission  Loss  Through 
Doubly-Periodic  3D  Elastic  Panels”,  Dey  S.  and  Shirron  J., 
Proceedings  of  IMECE  2006  ASME  2006  International 
Mechanical  Engineering  Congress  &  Exposition, 

(Nov.  2006) . 
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TECHNICAL  ISSUES 

•  Any  signal  processing  and  acoustic  beam  forming 
operations  strongly  depend  on  the  incident  wave 
front  remaining  relatively  undistorted  in  the 
presence  of  a  protective  window 

•  Guide  to  Perfect  Transparency:  Let  the  (solid 
material  Impedance)  (pcD)sol  =(pc)wat  (water 
impedance  ):  p=density  ;  c=sound  speed  ,  with  zero 
solid  shear  modulus  G 

•  ...NOT  POSSIBLE!  SOLID  SHEAR  MODULUS  “G”  IS 
ALWAYS  >0 

9  Therefore  different  kinds  of  structural  waves  can  be 
excited  during  the  fluid-solid  interaction 
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Governing  Equations: 


We  treat  a  general  class  of  problems  as  encountered  in 
structural  acoustics,  involving  a  submerged  elastic  plate 
that  is  subject  to  an  incident  harmonic  plane  wave  of  the 
form 


PinM^y^)=  Poe 


i  ( -  xk  v  -  yk  x-zkz+  tot ) 


Eq(1) 


Where  {kx  ky  kz  }are  the  acoustic  wave  number 
components  defining  the  direction  of  propagation 

n  =  i  kx  +  jky  +  kk. ) 

and  co  is  the  frequency  in  rad. /sec . 

The  thrust  of  this  paper  is  to  compute  the  acoustic 
transparency  lptr/pj  and  acoustic  reflection  /  prf/p0/ 


Governing  Equations  cont.: 

The  governing  equations  for  the  total  pressure  p  in  the 
acoustic  domain  and  displacement  vector  uj}  in  the  solid 
domain  (for  ei(ot  time  harmonic  response)  are  given  by: 


c  (d  p  !  dx  +  d  p  /  dv  +  d  p  /  dz  )  +  co  p  =  0 


and 


2  d  (oui.  2  d  /  j 

- ( - )  +  c  - ( - 

d  dx  j  dx-  $  dx-  dXj 


c 


where  c,  cd  and  cs  are  the  acoustic  wave  speed,  complex  solid 
dilatational  wave  speed  and  complex  solid  shear  wave  speed 
respectively. 

The  live  index  j  takes  on  values  j=1,2,3  and  repeated  indices 
sum  over  1  to  3. 


MATERIAL  CONSTANTS  WITH  LOSS  FACTORS: 

•  The  complex  wave  speeds  can  be  expressed  in  terms  of 
primed  real  wave  speeds  and  loss  factors  in  the  form: 

c]  =(c'd)2(l+irid) 
c2s  =(c's  )2(l  +  /»7?) 

Where  r)d  and  r)sare  the  dilatational  and  shear  loss  factors. 

•  Complex  wave  speeds  can  also  be  expressed  in  terms  of 
complex  E  modulus  and  complex  Poisson’s  ratio  v  as  : 

2  (l-v)£  2  _  0.5£ 

d  p?(l-2v)(l+v)  's  ps(  1+v) 

•  Upon  equating  the  corresponding  like  complex  wave  speeds 
we  can  easily  solve  for  complex  E  and  complex  v  in  terms  of 
c’d  ,c’s  ,  r|d  and  r\s  .  In  COMSOL,  the  isotropic  elastic  data  is 
entered  via  E  and  v  (in  MATLAB  syntax  format) , 

•  Thus  one  can  directly  enter  the  complex  E  and  v  in  the  form 

E=  E’  +  i*  E”  and  v=  v’  +  i*  v”. 


Governing  Equations  cont.  (scattered  formulation): 

•  COMSOL  has  a  feature  in  the  acoustics  module  that  permits 
working  with  the  scattered  (reflected)  pressure  formulation,  ps, 
rather  than  the  total  pressure  p  formulation 

•  However  the  formulation  used  here  still  uses  the  structure  of 
the  total  pressure  equations,  where  the  loading  of  the  incident 
field  now  enters  via  the  boundary  loading  at  the  incident  side 
fluid  structure  interface. 

•We  start  by  substituting  p  =  pinc+  ps  into  the  original  total  form 
time  harmonic  pressure  pde  to  get: 

c2(d2pJdx2+d2p/dy2+d2p/dz2)  +  w2P'  =  0 

U  LJ  l)  U 

•  Since  pjnc  is  a  solution  to  the  original  homogenous  pde,  it 
drops  out  after  the  substitution,  leaving  the  new  form  of  the  pde 
in  terms  of  ps  which  has  the  same  appearance  as  the  original  in 
terms  of  total  p. 


Governing  Equations  cont.  (boundary  conditions): 

•On  the  solid  interface,  (pjnc+  ps)=  the  total  pressure  is  applied, 
and  similarly  at  the  fluid  interface,  both  the  interface  plate 
acceleration  and  incident  fluid  acceleration  are  loaded. 

•On  incident  side  interface: 

inward  normal  acceleration  fluid  loading: 

nx_acpn  *u_  tt_acpn+ny_acpn  *v_  tt_acpn+p_i_acpr*i*(nx_acpn  *k 
x+ny_acpn*ky)/rhow 

normal  pressure  loading: 

-p*ny_acpn-p*nx_acpn-p_i_acpr*ny_acpn-p_i_acpr*nx_acpn 
On  transmitted  side  interface: 

inward  normal  acceleration  fluid  loading: 
nx_acpn  *u_  tt_acpn+ny_acpn  *v_  tt_acpn 
normal  pressure  loading: 

-p  *ny_acpn-p  *nx_acpn 

Where  p_i_acpr  =  Eq(1)  DRIVER  without  ei(ot;  rhow  =  po=1.0, 
kx=-kosin(0),  ky=kocos(0),  k=0  kQ=co/c . 


Governing  Equations  cont.  (radiation  boundary  condition): 

The  mesh  termination  at  the  end  of  the  fluid  domain  must 
include  some  sort  of  radiation  absorbing  boundary 
condition. 

a)  COMSOL’s  built  in  plane  wave  radiation  absorber 
(don ’t  need  acoustics  module) 

b)  COMSOL’s  built  in  PML  (Perfectly  Matched  Layer),  as 
indicated  by  the  extra  PML  zones  shown  in  Fig.1  above 
(do  peed  acoustics  module). 

c)  User  defined  Impedance  (Z=pjvn)  radiation  condition  , 
where  for  plane  waves  impinging  at  angle  0,  on  one  or 
more  completely  submerged  homogenous  elastic  layers 
and  using  Snell’s  law,  we  have  Z=  pfC/cos(G)  as  a 
radiation  boundary  condition  at  both  incident  side  mesh 
termination  and  transmitted  side  termination,  (don ’t  need 
acoustics  module). 


Governing  Equations  cont.  (fluid  periodic  boundary  conditions): 

•  Consider  next  the  boundary  condition  at  the  left  cut  vertical 
faces  and  right  cut  vertical  faces  (line  x=0  and  line  x=L  in  Fig.  1. 

•  The  response  at  the  left  cut  say ,  plcut ,  and  response  at  the 
right  cut  prcut,  is  not  known  in  advance,  however  we  do  know  a 
relationship  between  them,  namely: 

Prcut  =  Plcut^^ 

•Therefore  p,cut  and  prcut  are  not  independent  unknowns.  This  is 
a  “ Bloch-Floquet  theorem ”  type  condition  and  is  is  used  in 
refs. [5-7] . 

•  Applying  this  condition  in  COMSOL  is  not  intuitive.  There  are 
two  quirky  issues: 

a)  when  applying  the  periodic  boundary  involving  this  complex 
multiplier  condition,  the  user  applies  variable  p  on  the  left  cut 
and  p*exp(-ikxL)  on  the  right  cut  (the  conjugate  of  p*exp(+ikxL) ) 

b)  when  using  the  solve  parameters  advanced  settings,  one 
must  (according  to  COMSOL  staff),  check  the  box  engaging 
“Use  Hermitian  transpose  of  constraint  matrix). 


Governing  Equations  cont.  (solid  periodic  boundary  conditions): 

•  The  treatment  for  the  application  of  the  periodic  boundary 
condition  for  the  solid  follows  exactly  along  similar  lines, 
(where  upon  defining  displacement  components  as  u=u1  and 
v=u2)  we  have 

"rcut  =  uku/k*L  and  vra„  =  vlcu[eik*L 

Hi™  {  ^ lent  ’  V  lent  5  ^  rcut  ’  V  rent  ^ 

as  the  left  and  right  cut  displacement  component  values, 
analogous  to  the  pressure  left  and  right  cut  values. 

"The  displacements  are  enforced  similar  to  enforcing  p,  but 
instead  in  terms  of  COMSOL  variables  {u,v}. 

•Thus  {u,v}  are  applied  at  the  left  cut,  and  {u*exp(ikxL, 
v*exp(ikxL}  are  applied  at  the  right  cut. 


Post  Processing  Pressure  Fields 

•  For  the  scattered  formulation  employed  here,  on  the 
incident  side,  the  basic  COMSOL  variable  p_t_acpr  will 
actually  represent  the  scattered  component,  ps  ,  even 
though  it  is  labeled  as  total  in  the  post  processing  output 
list. 

•  The  actual  total  can  be  obtained  (if  needed)  by  simply 
adding  back  the  incident  pressure,  Eq(1),  (without  ei<ot). 

•in  the  case  of  the  presence  of  a  repeated  inclusion,  like 
in  the  Fig.1  model,  one  needs  a  strategy  for  computing 
the  transmitted  and  reflected  pressure. 

•  One  simple  measure  would  be  to  compute  integrated 
average  pressure  at  the  far  field  mesh  boundary  (but  not 
inside  the  PML  zone  if  present),  thus  getting: 
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Post  Processing  Pressure  Fields  (continued) 

"Next,  consider  an  alternate  method  to  process  a  representative 
transparency  pressure  in  a  variable  spatial  field,  namely  by 
computing  the  power  flow  across  the  y=constant  cut  boundary 

*  The  power  flow  can  be  computed  by  integrating  the  work  done 
over  one  time  cycle,  and  then  integrating  that  power/area  result 
over  the  top  (or  bottom  )  boundary  cut  surface  getting 


0 


•  After  normalizing  by  the  incident  wave  power,  can  be  used  to 
compute  the  “ pseudo  pressure  ratio”,  by  equating  the 
normalized  power  ,  to  the  normalized  power  in  a  pseudo  plane 
wave  (i.e.  ),  and  solving  for  the  transmitted  pseudo  pressure 
(and  similarly  the  reflected  value),  we  obtain: 


•The  power  relations  can  be  used  for  checking  the  consistency 
of  the  FEM  solutions.  Power  levels  can’t  be  created  greater  than 
the  input  normalized  incident  power  level  hence 


Scope  of  Applications: 

We  pass  through  a  sequence  of  increasingly  difficult 
problems  which  illustrate  the  application  of  the  periodic 
boundary  condition,  starting  with: 

(a)  a  simple  free  field  block  of  fluid  (no  plate) 

(b)  a  homogenous  submerged  elastic  plate  but  with  no 
periodic  inclusions 

(c)  a  submerged  plate  with  periodic  inclusions 
for  the  problem  type  shown  in  Fig.  1 . 


Free  Field  Propagation  Through  a  Simple  Fluid  Block 

*The  purpose  of  this  demonstration  is  to  illustrate  the 
enforcement  of  the  periodic  boundary  conditions,  with 
out  the  additional  complications  of  PML  absorbers,  and 
fluid  structure  interaction  boundarv  conditions. 


*  The  model  is  shown 
here,  where  the  bottom 
is  driven  with  Free  Field 
Plane  wave  pressure 
Eq(1 )  with  kz=0  and  ke=co/c, 
where  positive  kx  and  ky 
values  send  the  wave  in 
the  direction  shown 
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Free  Field  Propagation  Through  a  Simple  Fluid  Block  (Cont.) 

•The  solution  is  shown  here  and  perfectly  checks  out  against 
the  exact  solution. 


•  The  |p|  magnitude 
should  be  1.0  everywhere, 
and  the  direction  of  wave 
propagation  should  be 
along  straight  lines  that 
are  at  angles  0  =60°  to  the 
y  axis 

•  Where  it  is  noted 
that  these  wave  direction 
lines  are  perpendicular  to  the 
lines  of  constant  phase  wave 
fronts 


PHASE 


8  so 


-  0 


LOO 


LSO 


W-i  -J30 


■:  J.M 


MAG 


o  e 


o  o 


-  -  o; 


-  -  02 


0.1 


Sr-  0 


Plane  Wave  Propagation  Through  a  Simple  Flat  Plate 


"The  purpose  of  this  example  is  to  validate  acoustic-structure 
interaction  with  loss  with  the  enforcement  of  the  periodic 
boundary  conditions,  with  3  different  radiation  boundary 
conditions,  in  a  case  where  we  have  the  exact  solution. 

PR  ESSU  R  BPINC  A  POW  ERrPOW  INC  VS  A  NGLE  W  n"H  LOSS 


•  The  model  is  shown  below, 
except  the  inclusion  and  matrix 
material  of  the  plate  are  equal 
•This  is  basically  a  uniform  plate 
where  there  is  no  characteristic 
x  direction  length  ,  so  L=any 
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Damping  (Lossy  plate)  vs  No  Damping  Elasticity 
Solution  for  Submerged  2"  Steel  Plate 


Propagation  Through  a  Flat  Plate  With  Periodic  Inclusions 


•The  purpose  of  this  example  is  to  demonstrate  acoustic- 
structure  interaction  for  plates  with  periodic  inclusions, 
with  3  different  radiation  boundary  conditions,  for  the 


no  loss  model . 


•  The  FEM  model  is  shown  below 
except  the  inclusion  and  matrix 
plate  material  are  not  equal  j 
•This  is  non-homogenous  plate 
has  a  characteristic  period  in 
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Propagation  Through  a  Flat  Plate  With  Periodic  Inclusions 

•Here  we  show  the  amplitude  and  phase  of  the  spatial 
field  distributions  (no  loss  model)  in  the  front  and  back 
side  fluid 
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•  Lines  of  constant  phase 
illustrate  the  wave  fronts 
"The  magnitude  plots 
illustrate  the  expected 
Non-uniform  nature  of  the 

pressure  field  at  the  solid  .  .D  Q 

fluid  interface 
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Propagation  Through  a  Flat  Plate  With  Periodic  Inclusions 


"The  purpose  of  this  example  is  to  demonstrate  acoustic- 
structure  interaction  for  plates  with  periodic  inclusions, 
with  3  different  radiation  boundary  conditions,  for  the 


with  loss  model . 

•  The  FEM  model  is  shown  below 
where  the  inclusion  and  matrix 
material  of  the  plate  are  not  equal 

*  Note  that  the  sharp  peaks  are  „ 

now  smoothed  through  the  intro-  § 
duction  of  loss  (damping)  | 
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Propagation  Through  a  Flat  Plate  With  Periodic  Inclusions 


•The  purpose  of  this  example  is  to  demonstrate  acoustic- 
structure  interaction  for  plates  with  periodic  inclusions, 
for  a  variation  in  the  inclusion  stiffness  (with  loss) 
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•  3  different  inclusion  stiffnesses 
(i.e.  via  wave  speed  cD,  changes), 
are  compared 

oa 

*  The  stiffer  (bigger  cD)  the  inclusion, 
the  worse  the  transparency,  and 
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Propagation  Through  a  Flat  Plate  With  Periodic  Inclusions 

•The  purpose  of  this  example  is  to  demonstrate  acoustic- 
structure  interaction  for  plates  with  periodic  inclusions, 
for  a  variation  in  the  fluid  length  &post  processing  of  p 


•  2  different  fluid  model  lengths 
“Long”  (inset  below)  vs. 

“Short”  (i.e.  Long/10)  are  compared 

•  2  different  pressure  measures 
are  compared,  PAVG  vs  PPSUdo  y 

•  Short  model  PPSudo  consistent  |  0  6 
with  Long  model  PPSudo 
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THE  END 


